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SYNOPSIS 

RESIDUA! STRESSES IN PINE GRINDING 
S. BAIAIAH 

MASTER OE TECHNOLOGY 
Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 

In the present work an attempt has been made to evaluate 
the residual grinding stresses theoretically. Einite element method 
has been used to predict the distribution of residual stresses in 
the 'workpiece along the longitudinal as well as transverse direction 
during plunge-cut surface grinding. Eor simplicity only single pass 
grinding has been considered. Since the residual stresses are due to 
non-uniform thermal expansion of the workpiece, the temperature 
distribution in the workpiece during grinding has been evaluated 
following the model proposed by Sauer (9) with minor modifications. 

The theoretical results indicate that the peak residual 
tensile stress in the workpiece may reach the ultimate strength of 
the workpiece resulting in thermal failure of the workpiece. The 
results also indicate that grinding conditions have significant 
effect on the residual stress distribution. The peak residual tensile 
stress appears to decrease with decreasing wheel speed and increasing 
table speed, while increasing depth of cut increases its magnitude. 
The ' • - property of the grinding fluid also appears to have 

significant effect on the residual stress distribution. These 
results qualitatively agree with the published experimental 


results 



CHAPTER 1 


INTRODUCTION 


In grinding, the material Is removed by the action 
of many small cutting points of abrasive grains of uninden- 
tified geometry -which are held together in the grinding 
wheel by some bond material. Considering the action of the 
grain [Pig, (1)] it was shown (6) that at the beginning of the 
interference zone [Pig. (2)] , where a grain first contacts 
the workpieee surface, sliding and ploughing occur. At this 
location forces are not sufficiently large to remove the 
material. As the grain moves along the workpiece, it 
encounters a thicker layer of uncutmaterial , forces Increase 
and chip is removed. At the same time heat is generated at 
the grain-workpiece interface because of frictional sliding. 
Heat is also generated in the workpiece material because of 
plastic deformation. As the chip travels up the rake 
(leading) face of the grain heat is generated at the grain- 
chip interface. In adition to the temperature in the region 
of chip formation, the finished workpiece is also subjected 
to the wheel-workpiece interference zone temperature, as a 
result of this grinding burn or thermal cracks sometimes 
appear on the workpiece surface. 
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When any process which changes the shape or volume 
of an object in a non-uniform manner, the individual ele- 
ments of the object undergo different changes in their 
dimensions. Besides the forces in the workpiece which are 
in equilibrium with the applied forces, additional forces 
maintain the equilibrium within the workpiece or object 
after the release of external loadings. These internal 
forces cause stresses which are called "ilesidual stresses" 
or "Lo eked- in stresses". These stresses are maximum at or 
near the surface. Such a stress produces distortion, thus 
creating a compensating stress or resisting stress in other 
parts of the object. Since it is usually widely distributed, 
the compensating stress is much lower in magnitude than the 
primary residual stress and is therefore not significant. 
Once formed, the stresses remain unchanged indefinitely 
unless relieved' by plastic deformation or by failure. 

The basic causes of residual stresses are, 

(1) Structural transformation and consequent specific 
volume changes of micro-constituents in the workpiece. 

(2) Mechanical forces at the wheel -workpiece interface 
giving rise to inhomogeneous elastic-plastic deforma- 
tion of the workpiece. 

(3) Inhomogeneous thermal expansion of the workpiece 
resulting in thermal stresses. 

Under ordinary grinding conditions, the maximum 
temperature experienced by the workpiece will be well below 


the 
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structural allotropic transformation temperature of the wort- 
piece (9). So the possibility of residual stresses being 
generated due to allotropic transformation is negligible. 
Secondly, the investigation carried out by Mishra, etal.(25) 
reveals that the mechanical forces that are involved in 
surface grinding bear insignificant effect on the residual 
stress distribution. In the light of above facts it becomes 
evident that the residual stresses in grinding are nothing 
but the thermal stresses. 

-a typical subsurface temperature distribution in 
the workpiece for plunge-cut surface grinding is shown in 
Fig. (3). At any point below the surface the product 
(temperature x thermal coefficient of expansion x elastic 
modulus of the material) gives a good estimate of the 
resulting instantaneous compressive stress. Because of the 
colder underlying material restraining free expansion, the 
compressive stresses build-up as temperature ri ses. When 
the compressive yield, strength is exceeded plasticflow 
occurs causing upward flow on the free surface. After the 
material is returned to its average temperature the surface 
layer which had yielded is subjected to residual tensile 
stresses. The more the thermal upset, the higher is the 
magnitude of the residual stresses generated. 

It is well established that the effect of residual 
stresses have a similar effect on the fatigue behaviour of 
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materials as do mechanically imposed static stresses of the 
same magnitude. Thus significant residual stresses are 
beneficial if compressive and detrimental if tensile, parti- 
cularly in hardened materials. Thus for optimum design 
conditions, the control of residual stresses becomes an 
important factor in the control of fatigue strength of the 
component. Since the grinding operation is usually the 
last machining operation the analysis of residual grinding 
stresses becomes a relevant and a useful factor. 

1.1 PREVIOUS WORK: 

1.1.1 Grinding Temperatures : 

Temperature distribution in the region of interest is 
necessary for the calculation of thermal stresses. Several 
researchers in the past (1-5) have attempted to measure the 
temperature distribution in the workpiece during grinding. 
Because of the difficulties associated with grinding temper- 
ature measurement, attempts (6-9) have been made to evaluate 
the temperature distribution analytically. These temperature 
models employ the theory of moving heat sources proposed by 
Jaeger (10). The investigation carried out by DesRuisseaux 
and Zerkle (6) considers temperature in the viscinity of 
chip-formation and relates this temperature to the temperature 
experienced by the workpiece after grinding. The relation is 
established by considering the grinding geometry. Malkin and 
Anderson (7) divided the total grinding energy into chip 
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formation, ploughing, and sliding energy components. Their 
results indicate that virtually all the sliding and ploughing 
energies, and approximately 55$ of th e chip formation energy 
are conducted as heat into the workpiece. Both these models 
t have drawbacks in the sense that they do not take into account 

of the grinding fluid and the temperature rise due to the 
plastic-deformation of the workpiece material. 

The thermal model proposed by Sauer (9) on the other 
hand determines the workpiece temperature during grinding by 
combining experimental and analytical techniques. Some of 
the parameters have been evaluated by comparing theoretical 
and experimental results so that the best possible agreement 
between both temperatures is obtained. The cooling ability 
of the grinding fluid and the temperature rise due to plastic 
deformation are also evaluated. 

1.1.2 Residual Stresses: 

Several researchers (12-26) in the past have attempted 
to explain the origin of the subsurface residual stresses 
which result from the grinding of metals. The magnitude of 
these stresses have been measured by the experimental 
techniques such as parting out method, Layer removal method. 
Hole drilling method, X-ray diffraction method. Ultrasonic 
method (26) etc., in an attempt to reduce the incidence of 
fatigue cracks in engineering components. 
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Although there are numerous experimental results 
available on grinding stresses, very few attempts were made 
to evaluate these stresses theoretically (23-25). Nokano, 
etal . ( 24 ), have shown that workpiece is deformed into 
convex shape during grinding because of the temperature 
difference between the top and bottom surfaces, and that 
the profile of finished workpieces become concave. The 
residual stress distribution in the workpiece due to surface 
grinding has been obtained analytically by Mishra, etal. (25). 
In this work the temperature distribution in the workpiece 
was calculated for the case of strip heat source moving on 
the surface of a semi-infinite body. Finite element method 
was applied to evaluate the residual stress distribution 
following the initial strain method. No attempt was made 
to evaluate the effect of grinding variables on the residual 
stress distribution, or to establish the validity of their 
model with experiments. 

1 . 2 PRESENT WORK: 

In the present work an attempt has been made to 
evaluate the residual grinding stresses theoretically.lt 
has already been established that the residual stresses in 
the workpiece are due to the thermal stresses. In order to 
evaluate the temperature distribution in the workpiece during 
grinding, the model proposed by Sauer (9) has been adapted 
with minor changes. This model includesthe effect of grinding 
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fluid and plastic deformation and appears to yield results 
which are satisfactory. A computer program has been deve- 
loped to evaluate the temperature numerically. The data 
used were those given in reference (9) . 

Pinite element method has "been used to evaluate the 
residual grinding stresses along the longitudinal and trans- 
verse direction during plunge-cut surface grinding. Sub- 
parametric constant strain triangular elements (CST) were 
used. The analysis has been used to predict the effect of 
various grinding conditions on the residual stress distribu- 
tion. The results have been qualitatively compared with the 
experimental results published earlier. 



CHAPTER 2 


GRINDING TE5YH ERa IURES 


2.1 INTRODUCTION: 

In recent years several attempts have been made to 
analyse the mechanics of chip formation -with single abrasive 
grains. The model that appears to be most acceptable is the 
identation type model [Fig. 4(b)] proposed by Shaw,etal (2, 
30,31 ,32). They have obtained the average chip size and 
shape [Fig. 4(00 ] from the considerations of volume conti- 
nuity. It was further assumed that the grain tips are 
conical in shape which appears to be justified on the basis 
of available experimental results with single abrasive 
grains (31,32). Their equations for the Length of the 
chip, & c , the average width, b , and the maximum chip thickness, 
t,of an average chip are given by 


£ = /Dd 


h rt 


• V r .-r 
i w . d 

■ V — ~ v 


s 


D 


t 


4V 


w 


V r 


d 

D 


( 2 . 1 ) 

( 2 . 2 ) 

( 2 . 3 ) 


where D is the wheel diameter, d is the wheel depth of cut, 
is the workspeed, V g is the wheel speed, c is the number 
of active grains per unit area and r is the ratio of mean 
chip width to mean chip thickness. 
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'The variables c and r characterise the condition of 
the cutting surface of the wheel. Detailed discussion 
regarding these variables can be found in reference (34). 
Nakayama and Shaw -(30) have approximated the measured values 
of c by the equation, 

c = c Q (t - t Q ) (2.4) 

where c Q and t are constants for a given wheel. Sauer (9) 
has suggested this relationship in the following form 



b 1 


where t-j , m Q are constants for a given wheel. His evaluated 
values of these constants are given in Table 2 while the 
values of r are given in Table 3. 

The important factor in evaluating the grinding 
temperature is that a large amount of heat energy is generated 
deep in the workpiece due to plastic deformation that occurs 
beneath the grain. This heat source is called the "volume 
heat source". The detailed study of this aspect has been 
carried out by Sauer (9). In the present work his model 
and results have been extensively used for evaluating the 
temperature rise in the workpiece during grinding 

The strength of the volume heat source, q , can be 
assumed to be an exponential function of the depth below the 
surface, y, of the form (9) 
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■q v = q. e' 8y (2.6) 

V . r) 

■where q is the strength of the volume heat source at the 
v o 

surface (y = 0). The value of the decay coefficient, g , 
has to he determined experimentally. Sauer (9) has evaluated 
its value experimentally for AISI 1013 Cr steel and found it 
to he 27.73 1/in. 


In addition to this volume heat source there will he 
heat sources acting at the wheel -work interface and at the 
chip-work interface , both of these sources are called "surface 
heat sources". The plane over which these heat sources 
extend is shown hy the line B T C' in fig. t 4<b)3 . Its 
width, hg., equals the width of the chip, b , and the length 
is represented hy and is equal to B’C' . Thus, 


h 


G 



(2.7) 


The approximate relation for is given hy (9) 


SL 


G 


3.33 


i 

i 

i 

v 



M_ TT H r 
c w 


11 

t 


(2.3) 


where 3? is the horizontal component of grinding force, M Q 
is the total number of grains in contact with the workpiece 
at any time , is the warm hardness of the workpiece 

material and t’ is the instantaneous chip thickness. Approxi- 
mating the chip shape to he a slender wedge (6) instantaneous 
chip thickness, t’ , can he evaluated from 
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(x + £ ) 

tt = t (2.9) 

2 £ c 

where x is the horizontal location of the maximum chip 
thickness from the origin 0. [ Pigs. (5 ,6 ) ] . 

2.2 ANALYSIS: 

The workpiece can be represented mathematically by 
a region R as shown in Fig. (5). 

| -’.co < X < co j 

i 

1 

R = *-co< Z < co ' (2.10) 

j 0 <- y < co 1 

The area of contact between the workpiece and the wheel is 
from x = — £ to x = £ and infinitely far in the z direction. 

Thus we can assume that there is no temperature variation in 
the z direction when the workpiece moves with the constant 
velocity V relative to the stationary wheel in negative x 

W 

direction. The coordinate system is fixed to the wheel 
having 0 as the origin as shown in Fig. (6). Since the depth 
of cut, d, is much smaller than the length of contact, the 
plane A300D can be approximated as y = 0 plane. 

The model of a semi-infinite body can be used to 
simulate the finite body if the body is large enough so 
that end effects can be neglected at the centre portion of 
the body and the actual teraperature distribution can be 
approximated by the steady state temperature distribution. 
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2.2.1 Temperature Field in a Serai-infinite Solid due to 
Moving Sources of Heat : 

The equations for calculating the temperature distri- 
bution in the workpiece during surface grinding have been 
simplified by Sauer (9). In deriving these equations, the 
theory of moving heat sources for the case of line heat 
source proposed by Jaeger (11) has been used. 

The* temperature distribution in the semi-infinite 
body due to moving line heat source of strength q^ can be 
evaluated from (11) as 


^ L (x,y) 


2 -(y-y* )‘ 

r x - ( *-x'+v t) e *p o 4 ' et - ■ 


= [ — J exp ( 

pc 0 


4 k T 


} I 


■} 


+ 


exp { 


-(y+y’ )' 

4 k x 


•} 


4 7tk T 


] dx]~ [ 


f 


4tt k t 
h/K 


2 Pc o /if k x 


y+y’ +2 § k x 

erfc{ } exp{ (y+y') | + k(|) ? 


2 /kt 
x’ + v’x)‘ 


4 kx 


dx ] 


( 2 . 11 ) 


q^ dz' dx 1 


where x’ , y', z' are the coordinates of the point source 
of strength q p acting at time t t f P i s the extent of th 
plastic zone as shown in Pig. (4(a)), v' is the velocity of 
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the heat source, k is the thermal diffusivity of the work- 
piece material, £ is the thermal conductivity of the work 
piece material and h is the convective coefficient of the 
cooling medium. 

Sauer (9) has simplified the above equation by 
introducing the following non-dimensional quantities, 


X 


2k 


v 

y = y 

u = 


H = 


<S> = 


L = 


k 2k 

K 7 

Tt K v' 


2q, 


v 


T 


Z 

2k T 


( 2 . 12 ) 


where X is the dimensionless x coordinate, X is the dimension- 
less y coordinate, U is the dimensionless time variable, H is 
the dimensionless convective coefficient, 4> is the dimension- 
less temperature, L is the dimensionless half length of the 
heat source, T is the temperature variable and q Q is the 
source strength parameter. Equation (2.11) thus takes the 


form, 
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-(X-X'+U)' 


-( Y-Y’ )' 


~ J exp { } [ exp { } 


% 8k o 


2U 


2U 


+ exp {■ 


-(Y+Y t ) 
2U 


} 1 


dU 


v H 


U\ q o 4k 


7T 

2 


00 Y+Y’+HU -(X-X 1 +U) „2 

J erf c { — -} exp { } exp { H(Y+Y* )+ 


/'2U 


2U 


/ T J 


(2.13) 


This equation can he considered in two parts, the first part 
representing the temperature rise due to the heat source 
without convective loss from the surface and the second part 
representing the temperature reduction due to convective 
cooling from the surface. 


2.2.1(a) Temperature Rise due to Single drain Surface 
Heat Source: 

Assuming the strength of the heat source, q^, 
'to be uniform over its length, wt -c* 


= q 0 dx’ for - -<.x’ <l $, ; y’= 0 ( 2 . 14 ) 

Integrating the first part of equation (2.13) over the extent 
of the heat source, th e temperature rise, 4> ^ , will be 


$ 


= h 


JT 

2 


exp(- 


str 


/u 


{ erf c( ^tS) - erfc(M±2) 


dU 


/~2U 


»2U 


(2.15) 
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2.2.1(b) Temperature Rise due to Single G-rain Volume 
Heat Soux-ce: 

The strength of the volume heat source can be 
evaluated from (9) as 

q^ = 8 e ^ ^ dy’ dx ! for - £<^ x’ < H 

0 <_ y' < co (2.16) 

where 8 is the decay coefficient. Introducing the dimension 
less decay coefficient, n , as defined by 



and after integrating the first part of equation (2.13) over 
the extent' of heat source the temperature rise due to volume 
heat source, is given by 

00 

I l erfe( «±2 ) . e rf=( 2±S±E ) } 

8 o /2U /2U 

exp( n 2 5 >{ exp( qY) erf c( ; 

d /'2U 

+ exp(-nY) erfc( } dU (2.17) 

/ 2U 

2.2.1(c) Temperature Rise due to Interference Zone 
Surface Heat Source: 

Assuming the surface heat source, q^, to be lin- 
early varying with x T , we get 
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q^ = q Q (1 + y) dx’ for - 2, < x' < £ 

y ! = 0 ( 2 .1 S) 

After integrating the first part of equation (2.13) over the 
extent of heat source, the temperature rise due to interference 
zone surface heat source, <J > ^ , is 

4>5 = yr- I ex P(" lu> [ 4 -gjCx+I+U) { erfc(^~^) 

O 

- erfc(S±£±D ) > + { exp( AM±uf)-exp(-(i|±ial))l a u 

/ 2TJ 

(2.19) 

Temperature reduction due to surface cooling: 

After integrating the second part of equation (2.13) 
over the extent of the heat source, the reduction in tempera- 
ture, <J>c» due to surface cooling is given hy 

00 2 

A = f exp (HI + 2_2) erfc(^ip2) [ (X+L+U) { erfc( - ~ - -r~) 

^ 4 1 o ^ /2U /2U 


.erfo^)} + /f { exp(A ^ al!)-exp(. -iI i §± a)! } ] an 

( 2 . 20 ) 

Hence the effective temperature rise due to interference zone 
surface heat source is given hy 



4> 




3 


5 
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2.2.1(d) Temperature Rise due to Interference Zone 
Volume Heat Source: 

The source strength is assumed to he linearly 
varying with x T so that 

q L =gq (1 + — ) e" Sj ’ dx’ dy’ for - £< x' < ~ 

0 <. y ! < 00 

( 2 . 21 ) 

After integrating the first part of equation (2.13) over 
the extent of the heat source, the temperature rise, 4> 4> due 
to interference zone volume heat source is 

$ 4 = J exp(— ^p2) [exp( nY) erfc( i -~ ± -^) 

4 SL o d / 2U 

+ exp(- ny) erf c(-^— L—-^) 1 t (Z + L + U) 

*2U 

{erf0 (M± 2) . erfo(Ai±U) } +/ '£ { SIp (lM £) 

/2U /2TJ 77 2U 

2 

_exp( ~C^r^±Sl- )} ] du (2.22) 

Temperature reduction due to surface cooling: 

After integrating the second part of equation (2.13) 
over the extent of the heat source, the reduction in temper- 
ature, <|>g, due to surface cooling is given hy 



IS 


*6 = 1 H / [ ( T-t-L+U ) (erfc(^) - erfc(^) } 


✓ 2U 


/2~U 


+ 


/f { e X p(li ™!) - expCdlzW) } 


{ 


exp(HY+ — ^) 

__ 


e rfc (2±32) + .“>< " Y+ ^ 


/ 2U 


H - n 


erfc(^ r--- -) > ] dU 
/ 2U 


(2.23) 


Thus, the effective temperature rise due to interference zone 
volume heat source will "be 

^ 4,6 = ^4 “ ^6 

and the overall temperature rise of the workpiece is 


*-*1 +,f 2 + V,5 + 4, 4,6 


To obtain the dimensional temperature, T, the nondimensional 
temperature,^ , has to be multiplied by the factor 
2q k 

( o_ ) 

11 K v w 

The characteristic parameter, a, for different heat sources 


are (9) 

Vi) 


“ 1 T 

T-yvr (1 + f> 

G- Or C 


l o<2) 


“2 Q P , X 

(1 + — ) 


”8 M c 
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• 

01 1 % 

q o(3,5) £ c B 

a 2 

q °(4,6) £ 5 

° (2.24) 

where aj is the fraction of the total grinding energy 
entering the workpiece through the surface (y = 0),a 2 i g 
the fraction of the total grinding energy generated in the 
workpiece due to plasti c deformation , Q-, is the rate of 

total grinding energy generated, and 3 is the work width. 

a = + a ^ and Qp is given hy 

F V 

Q = — ^ Btu/sec. (2.25) 

■“ ( 1 2) (77S) 

The temperature so obtained is the steady state tanperaturc 
over the ambient temperature. 

In deriving the above equations the thermodynamic 
properties of the material of the semi-infinite region 
are assumed to be independent of the temperature variation. 
But for the real materials it is not true. However, 
comparison of temperature distribution with constant 
properties and with variable properties indicates no 
significant difference (9). Hence the average value of 
the material properties have been used in the present work. 
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As per the theories of hardness test (32), the warm hard- 
ness of the workpiece material has teen taken a.s three times 
the yield stress of the workpiece material at its average 
temperature during grinding. 

Using equations (2.13) through (2.22), the temperature 
rise can be evaluated. The details of the numerical technique 
used for evaluating this on computer is discussed later. 


CHAPTER 3 


RESIDUAL GRINDING STRESSES 


3.1 INTRODUCTION: 

As discussed earlier the residual stresses in grinding 
are the thermal stresses generated due to non-unif orm thermal 
expansion of the -workpiece. In the present work, the average 
values of the thermodynamic and mechanical properties of the 
workpiece material have been used and the material is assumed 
to be isotropic. The strain-stress relations for a three 
dimensional thermal stress problem of an isotropic solid are 
(36), 


e x - a t T = I ta x " v(a y + cr z )] 

£ y " a t T “ I [a y “ V(a x + a z )] 

e z " a t T = h 0 z " v ( a x + a y )3 

Y = I £L . Y = Tyz • v = T xz 

xy G ’ yz G * T xz G 


(3.1(a)) 


(3.1(b)) 


where a . a rr . a _ are the normal stresses along x, y, z axes 
x y z 

of the coordinate system [Pig. (5)3 , a 3 are the normal 

strains along x,y,z axes, t ,t , are the shear stresses, 

Y „ , Y , y are the shear strains and v E, G, n are the 
poisson - ratio, elastic modulus, shear modulus, thermal 
coefficient of linear expansion of the workpiece material 
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respectively. Equations ( 3. 1 (b) ) are not affected by the 
temperature as free thermal expansion does not produce 
angular distortion in an isotropic material. The volume 
expansion, e, and the sum of the normal stresses, 0 , are 
given by 


e _ e + e + e 
e ~ x + y z 


(3.2) 


0 = cr x + cr y + 0 Z 


e and e are related by the equation, 


(1 - 2V) 
E 


(3.3) 


Using equations (3.1(a)), (3.2) and (3.3), we get 


e = | (1 - 2v)0 + 3a t T 


(3.4) 


Tnus, the stresses are 


a = A e + 2G e 
x x 


a E T 
1-2V 


a . E T 

cf ss A g + 2G £ - — — — — 

y y i -2 v 


a = A e + 2G£ „ 
z z 


a t E T 
1 -2 v 


(3.5) 


where 


A = 


v E 


(1 + v)(l_2v) 

E 


(3.6) 


G = 


2(1 +v ) 
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The equations of equilibrium for an elastic continuum without 


body forces are given by (36) 


3b 


x 


3x 


3 x 


j£L 


3t 


xz 


a y 


a z 


3b 3x 3 t 

L + SL + I* 

ay a x a z 


3b 3 t 3t 

z xz yz 

a z + ax + ay 


o 

o 

o 


(3.7) 


and the strain-displacement relations are 

e a u . e _ a v . _ a w 

x 3 "x 5 y 3 y ; £ z ” 3 z 



3u 

DV 

xy 

3 y 

+ a x 


a u 

a w 

xz ” 

Oz 

+ a x 


a v 

a w 

yz 

a z 

+ Ty 


(3.8) 


where u,v,w are the displacement components along x, y, z 
axes respectively. Using equations (3.1), (3.7) and (3.8), 
the equations of equilibrium can be expressed in terms of 
displacements as 


( ^ + G) -g— : + G 7 2 u 

( X + G) + GV 2 v 

( X + G) 4| + G v 2 w 


a E 
t 

a i 

1-2 v 

a x 

“t E 

a t 

1-2v 

a y 

a t B 

a t 

1 -2 v 

a z 


= 0 


= 0 


= 0 


(3.9) 
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where V = 


8 4 - 3 

TIE + 


3 7 


3 z 


Any particular solution of equation (3.9) will reduce 
the thermal stress problem to an ordinary problem of surface 
forces, l'he solution for u,v,w, subject to the boundary 
conditions of a particular case will give, by means of 
equations (3.S) and (3.1), the values of the stress components 
of the corresponding problem. One way of finding particular 
solution to equation (3.9) is to take 


u 


JiL 

3 X 


3y ’ 


w 


3w 
3 z 


(3.10) 


where the thermoelastic displacement potential, ip , is a 
function of x,y,z and also a function of time if temperature 
varies with time. Using equation (3.6), equations (3.9) 
can be rewritten in the form 

3 e + (i.2v ) V 2 u = 2(1 + v ) a, 3 T 


T5E 


t 3 x 


Ty + (1 ~ 2 V) 2(1 + V > a tff 

Af + (1 - 2v )V 2 w == 2(1 + v) 


(3.11) 


0 . „ 3 u 3 v 3 w 

Since © = + 3 ~ + yi » 

2 

e = V if; 

and equation (3.11) becomes 


equation (3.10) will yield 


(1 -v)- 


3x 


(V ip )_ 


(1 + V) a 


3 1 
t T*x 
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0 - V) k (V *> = (1 + V) a t f-f 


ay 


(i -v^CvtJ.d +v)OH! |_S 


t 9 z 


( 3 . 12 ) 


Equations (3.12) arc satisfied by taking the function 4> as 
a solution of the poisson’ s equation, 

7% = ( ) “ t T (3-13) 

Now any thermal stress problem can be solved by solving the 
above equation provided the temperature function; T(x,y,z), 
is known. 


3.2 ANALYSIS: 

The temperature distribution in tne workpiece during 
grinding, as discussed earlier, is independent of the 
z direction. Thus equation (3.13) can be further simplified 
to the plane— strain condition. 


2 ? 

9 y 3^ 

+ 

9 x 9 j 


= ( T - $ a t 


(3.14) 


Since the temperature distribution is evaluated numerically, 
numerical solution method has been used for evaluating thermal 
stresses also. The finite element method has been used for 
evaluating the thermal stresses since it yields better results 
than the finite difference method. 


Solution of Poisson’s equation by Finite element method : 
The equation (3.14) can be written as 


if. 

9 x 


4 

T 


+ 


9 2 4> 
9 y 


+ o 


= 0 


(3.15) 



26 


where 0 = -( ^ I (3.16) 

Equation (3.15) is the necessary condition for functional 
E to he minimum. This functional E, is given hy (38) 


E 


//[*( 

R 


^L) 

a x' 


2 



Grp ] dx dy 


(3.17) 


where R is the region of considers tion^Eig. (5 )] - Thus the 
solution, ip , of equation (3.15) is that value of ip "hich 
minimises the functional E [equation (3.17) 3 . Minimisation 
of the functional, E, and the equations for the element 
stiffness matrix and the load vector are clearly given in 
reference (38). Einal form of those equations are 



, , a n. a n. a e, a R. 

/ J + ~ y g y' ) ( ^ rl 

A e ^ 


(3. IS) 



f G N. dA 

A (.) 1 


(3.19) 


( 6 ) ( 6 ) 

where k. . and p. are the elements corresponding to ith 

row and jth colu an of the element stiffness matrix, (k^ \ 

(o') 

and element load vector, {p }, respectively. , E. are 

( e) 

shape functionsassociated with the nodes and A K ' is the 
element area. The assembled equation is 


[£] W + { P> = 0 


( 3 . 20 ) 


where [ K] is the assembled stiffness matrix, {>} is the 
solution vector of the nodal values of , and (E) is 
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the assembled load vector. The set of simultaneous 
equations (3.20) are solved by Gauss elimination method 
to obtain the solution vector, {$} , subject to the 
boundary conditions discussed below. 

B oundary conditions : 

Referring to Rig. (7), the workpiece in an actual 
surface grinding is fixed on the magnetic chuck. The experi- 
mental results that are available in the literature (21) 
reveal that the depth upto which the residual stresses are 
dominant is very low compared to the total depth of the 
workpiece. Hence, at depth, say y = y^, in the material 
itself the displacements are very low. By trial and error 
method the depth y^ was found to be 0.14 in. , where tne 
displacements are insignificant compared to the displacements 
near the surface. This is achieved by prescribing the value 
of 'P to be zero at y = y^. The displacements are computed 
and compared with the displacements in the upper region, 
i.e. , y< y^. It was found that the displacements at 
y = 0.14 in. were of the order of 10 to 10 whereas the 
displacements near the surface y = 0, were of the- order of 
10-1 

The continuum is divided into subparametric triangular 
elements. A linear interpolation function is used for the 
description of the geometry and a quadratic interpolation 
function is used for the variation of ip within the element. 
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A single degree of freedom system is sufficient to satisfy 
the boundary condition- The temperature distribution within 
the continuum is specified by giving as input the temperatures 
at each of the nodal points. The variation of ip within the 
element, ip has been taken as 


ip 


U) -t * 1 { * nodal } 


( 3 . 21 ) 


Displacements : 

The displacements are obtained from 


u 


v 


< 3 x 

< 

■ li_ 
! 3 y 


( 3 . 22 ) 


where 


3ip 

"5x 


! 


D 2 3N 3 3N 4 3% 


3 x 3 x 3 x ax 


3 x 


3N 


6 * 


3 x 


*1 
ip 2 

*4 

*5 
* £ 


dip 

3y 


| 3lD| 

Nr 


a No 3 % 3 F 


a y 


3 y ay 


3 IT C 

TF 


3N e 

TF 


tp 

b 

^p 

Ip 

Ip 

Ip 


4 


6 ■ 
(3.23) 
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St rains : 

j x .11 — • 

The strains are given by 
. 2 


'X 


3 i|> 

T“2~ 

3 x 


xj 


2 

3 tb 

, 3 2 t 
' 3x 3y 


2 

3 il> 


(3.24) 


r» 


2 , 


are 


Tor plane strain condition where 

3 x“~ 3 y*" 

obtained by differentiating 9^/ ®x, 3 tp / 3y with respect 
to x and y, respectively in the same manner as 3v/^x, H'/3y 
are obtained from equation (3.21). 


Stresses : 

Trom the ' strain-stress relations given by equationC 3. 1 ) , 
for the case of plane strain condition, the stress components 
can be expressed by 


{a} - 



-iDl \ {e> - {e ) 
* U 
\ 


where {e} = 


x 


\ e y 


{ Y V j 


and ( e Q > = ( 1 + v ) 


(3. 25) 


I « t 

t 




T 


T 


(3.26) 



and 
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j 

v/(1-v) 0 j 


0 (1-2 v) ■ 

These stress components can also be obtained from equation 
(3.5). a , the stress along the transverse direction is 
obtained from equation (3.5). The strains and the stresses 
are constant within the element. 

The element stiffness matrix and element load vector 
are derived for the subparametric triangular element. These 
matrices along with the ^ -displacement and displacement-strain, 
transformation matrices are given in Appendix. Boundary con- 
ditions are incorporated by eliminating the rows and columns 
of the assembled stiffness matrix and the rows of the assembled 
load vector corresponding to the zero degrees of freedom. 

Now having obtained the temperature distribution and 
using equation (3.14) thermal stresses can be evaluated 
numerically. The details of numerical method used for evalua- 
ting this on computer is discussed later. 


EDI = 


E(1 - v) 
(l+v)(1-2 v) 


1 

v/(1- V ) 
0 



CHAPTER 4 


NUMERICAL SOLUTION 


The computer program is divided in two parts. The 
first part calculates the grinding temperatures numerically. 
The temperature data so obtained are used as input to the 
second part to evaluate the thermal stresses. 

Equations (2.15) through (2.23), which describe the 
steady state temperature distribution in the moving semi- 
infinite solid, cannot be integrated analytically . Hence 
to obtain numerical results for a specific problem, a 
general purpose computer program has been developed. The 
numerical integration is carried out by using 10 point 
Gaussian-quadrature. A standard subroutine is used for this 

purpose. The results so obtained are compared with the 

results given by Sauer (9). Good agreement is evidence'’ 

from Eigs. (8 through 13). Since the computer is limited in 

size and the numbers which it can handle, two problems that 
arise in this solution method had to be overcome. 

Some of the equations contain terms of the form 

e (s -b ) erfc(s) (4.1) 

Eor large values of s (say > 9) the exponential part of 
the expression (4.1) increases rapidly beyond the limit of 
computer (10 for IBM 7044), while the value of the complete 
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term remains 
complementary 

erfc(s) = e“ s 


finite for all possible values of s. Hence, 
error function had to he approximated as 


5 

I 

i=1 (1+p s) 1 


a i 


for s > 0 


(4.2) 


the 


This is based on the formula given by Hastings (39). Here a i 
and p are constants and have the values given below. 


p = 0.3275911 
a 1 = 0.22583SS4 
a 2 = -0.25212866 
a 3 = 1.25969513 
a 4 = -1.28782245 
a 5 = 0.94064607 


Using equations (4.2) and (4.3), expression (4.1) can be 
calculated by defining a new function, 

2 5a. 

experf(s) = e s erfc(s) = I i — T for s>0 (4.4) 

i=1 (1+ps) 1 


Thus, 

( 2 -h 2 ^ -b 2 

e^ s ' erfc(s) = e~ experf(s) (4.5) 

where s and b are related such that 


s 2 < b 2 for s < 0 (4.6) 

The new function is not defined for negative value of S but 
this is not necessary. Hence both parts of expression (4.1) 
remain constant for all finite values of s. 
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The second problem faced during computation was 
<1 

due to the definition of finite upper li'iit for the inte- 
gration of equations (2.15) through (2.23), The upper limit 
00 implies that the calculated temperature is steady state 
temperature. However, this amounts to saying that the 
upper limit of integration is the time, say U , at which 
the point of the coordinate system has reached its 93 , 99 % 
of the maximum temperature. This time, IT , is a function 
of the half length of the heat source, L, and the position 
of the point in the coordinate system. The expression for 
U ss as given by Sauer (9) has been used in the present 
program. 

IT = 2L - X + 0.5Y + 10 (4.7) 

o o 

Thus the steady state temperature distribution has been 
obtained by using U as the upper limit of integration. 

For time less than U__, transient temperature field can be 

S o 

calculated. 

The second part of the program calculates the thermal 
stresses numerically using equation (3.14) by employing the 
finite element method. For checking the correct assemblage 
of the assembled stiffness matrix and the total load vector 
etc., a sample problem with 4 elements was solved and the 
results obtained were checked with hand calculations. 

The flowcharts for both part one and part two of the 
computer program are given in Figs. (IT. 1) and (¥.2). 


START 


/' READ A1TD PRINT V w , Vg ,d ,D,H w ,k ,K, c* 2 , ' 

B, F , h , GRAIN SIZE, WEED GRADE AND NODAL 

P COORDINATES 

_ — - 


CHOOSE r, t 1 ,m , CORRESPONDING TO GRAIN SIZE 


AND WHEEL GRADE 


DO 1=1 , NUMBER OF NODES 


V/ 


calculate non-dimensional i aR-aiems x , y , l , h / , $ ; a 

AND DIMENSIONAL PaR/AIETIRS t,c,M p 


IS X<(-L)> 


\T NO 


YES 


=0 > 

■ — - ! 


CALCULATE V^G’S 


/IS 1 X I > L> 
'• ^ 

NO 


YES 




OAlOUUlTE Vl) (I) ’V2) (I) ’V3,5) (I) ’ q o(4,6) (I) 


COMPUTE ^ ,^2^3»^4^5^6’S,5 ^4,6^ 3Y USITTG 
10 POINT GAUSSIAN QUADRATURE AND CALCULATE 
DIMENSIONAL TOTAL TFMPERaTU IE T 



( PRINT AND PUNCH NODAL 
TEMPERATURES 


V 


/ 


( STOP 


FIO 1 : Computer flow chart for caxculation of 

temperature distribution In the worhpieoe. 



/ READ itND PRINT, NODAL COORDINATES (x,y), \ 
i GLOBaL NODE NUMBERS CORRESPONDING TO \ 
! ELEMENT NODES, NODaL TEMPERATURES , j 
\ NUMBER OE ELEMENTS, <* t , v, E. ! 


initialise the Matrices [k], ip 3 j 


! DO KK = 1 , NUMBER OP ELEMENTS (NE) j 


1 PIND OUT, THE GLOBAL NODE NUMBERS j 
I CORRESPONDING TO ELEMENT NODES AND I 
! THEIR COORDINATES ; 


COMPUTE t k^l , { > S 


! ASSEMBLE THE MATRICES [El, E P 1 AND 
| INCORPORATE THE BOUNDARY CONDITION 


P SOLVE THE SIMULTANEOUS EQUATIONS 

! K {ip} = ( P } TO OBTAIN THE SOLUTION 

I VECTOR W 


COMPUTE NODaL DISPLACEMENTS, ELEMENT 
STRAINS, STRESSES AND VERTICAL DISTANCE, 
CTRY, OE THE ELEMENTS CENTROID BELOW THE 
ORIGIN 

J 

■ — r 

i 

PRINT CTRY, a x » RESIDUAL STRESS ALONG ' v 

LONGITUDINAL DIRECTION, a z , RESIDUAL j 

STRESS ALONG TRANSVERSE DIRECTION ! 

CORRESPONDING TO THE ELEMENT KK. 


/IS 

/• KK = NE 
< 


YES 


S' 

^ NO 


STOP \ 


PIG.(N.2): Computer flow chart for evaluating tEe 
residual stress distribution by finite 
element method. 



CHAPTER 5 


RESULTS AND DISCUSSION 


As discussed earlier, the semi-infinite model has been 
used for calculating the steady state temperature distribution 
in the finite workpiece. Referring to Eig.('U)* the steady 
state temperature distribution so obtained decays very fast 
with the depth below the surface . It becomes almost negligible 
at a depth approximately equal to 0.075 in. ( =0.2 mm). This 
is sufficiently small compared to the total depth of the 
workpiece to justify the simulation of the finite workpiece as 
the semi-infinite body. 

A typical residual stress distribution in the work- 
piece obtained by using the finite element method is shown in 
Pig. (15). It indicates that the layers very close to the 
surface are under high tensile stress. The stresses vary 
from tensile to compressive and the magnitude of these stresses 
decrease with increasing depth below the surface. However, from 
the experimental results available in the literature (14-21) 
it can be said that the nature and distribution of the residual 
stresses may vary widely with the grinding conditions. Positive 
values indicate tensile stresses and negative values indicate 
compressive stresses. Investigators in th<= past (14-21) have 
reported two stress values, one along the longitudinal direction 
and another along the transverse direction. It was shown that 



35 


the stresses in other directions will have intermediate 
values that lie between the two stress values mentioned. 

Hence, the distribution of residual stresses in the longi- 
tudinal and transverse direction are only considered and are 
shown in Pig. (15). A typical residual stress distribution 
measured by Letner (21) is shown in Pig. (16) and a schematic 
residual stress distribution for the case- of dominant thermal 
source as given by Colwell, etal* (19) is shown in Pig. (17). 

A good qualitative agreement can be seen between the se figures. 
If the peak residual tensile stresses which occur at or near 
the surface reach the ultimate strength of the workpiece 
material, the thermal failure in the form of "grinding cracks" 
can occur. Thus, it can be concluded that the magnitude of the 
residual tensile stresses have to be minimised to avoid thermal 
failure of the workpiece. 

The grinding parameters such as table speed, wheel 
speed , cooling ability of grinding fluid, depth of cut, hard- 
ness of workpiece material, grain size, wheel grade etc. have 
a direct influence on the residual stresses generated during 
grinding. Effect of these parameters on residual stress dis- 
tribution had been studied in the past experimentally (21) 
to make the grinding operation more useful. This has been 
evaluated theoretically in the present work. 

Referring to Pig. (IS), the magnitude of peak residual 
tensile stresses decrease with increasing table speed. This 
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is "because the fraction, a , of the total grinding energy 
entering the workpiece decreases (Table 1) as the table 
speed increases. Thus, the maximum tenq erature to which 
the workpiece is subjected decreases. On the other hand, 
at low table speeds, the reduction in the magnitude of 
residual stresses appear to be insignificant. At low 
table speeds the cooling property of the grinding fluid 
becomes so Qffective that it overcomes the effect of reduced 
table speeds. As a result of this, the temperature to which 
the workpiece is subjected is low. The obtained variation 
of residual stresses with table speed [ Pig. 13(a)] can be 
compared with the variation measured by Plusch (21) 

[Pig. 13(b)! which shows good qualitative agreement. 

e^ppeoa/s "to (UiCcXfe 

Table 1 that the fraction, a , decreases 

with decreasing wheel speed. Thus the magnitude of residual 
stresses increases with increasing wheel speed. This is 
shown in Pig. (19). Pi eld and Kahles (33) also found the 
residual stresses to increase with wheel speed [ Pig. 19(b)]. 

As the depth of cut is increased the fraction, a , 
increases. Consequently the maximum temperature at a point 
in the workpiece increases. Thus the peak residual tensile 
stresses increase with the increase in depth of cut. The 
effect of depth of cut on residual stress distribution as 
shown in Pig. (20) indicates this to be the case. One way of 
measuring the residual stress is by removing successive layers 
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from the work-piece surface and measuring the change in the 
curvature resulting from release of the stresses in each 
layer (12). letner (12) has given the relationship between 
the magnitude of residual stresses and the corresponding 
curvature changes. In the experiments carried out by £ubsh(21), 
curvature change was found to increase with increasing depth 
of cut [Pig. (21)]. The re-sults of Pig. (20) and Pig. (21) seem 
to match qualitatively. 

It can be seen in Pig. (22) that the magnitude of 
peak residual tensile stress decreases with increasing 
convective coefficient of the grinding fluid. After reaching 
a certain value of convective coefficient this reduction in 
the magnitude of residual tensile stress is insignificant. 

This is because, the reduction in magnitude of peak residual 

conv€e"tcv€ 

tensile stress depends upon good propc rty of the 

grinding fluid, while large va.riations in the cooling property 

njR>o ooo5 bUl/'c UA see °f)I . , 

Thas only a very small effect on the temperature rise and hence 

on residual stresses. 

Malkin and Anderson (35) have shown that as the 
grain size (nominal grain diameter) is decreased, the- number 
of active grains per unit area increases. Therefore, the 
strength of the single grain heat source must decrease with 
decreasing grain size (increasing mesh number) resulting in 
lower heat energy generation which affects the residual 
stress distribution. Results shown in Fig. (23) also indicate 
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that the residual stress values decrease with decreasing 
grain size (increasing mesh number). This aspect, nowever, 
needs further investigation. 



CHAPTER 6 


CONCLUSION' 


The theoretical model using the finite element 
technique for evaluating the residual stress distribution 
in the workpiece during grinding appears to yield results 
which agree qualitatively with the published experimental 
results. It is observed that the magnitude of the peak 
residual tensile stress can reach the ultimate strength of 
the workpiece resulting in thermal failure of the workpiece. 
The grinding conditions have significant effect on the 
residual stress distribution. 

The peak residual tensile s cress appears to decrease 
with increasing table speed, decrease with decreasing wheel 
speed and increase with increasing depth of cut. The 

to nve<£Eve 

; property of the grinding fluid also effects the 
residual stress distribution significantly. 

G * ! Future Research 

In the present work a theoretical model has been 
presented and it needs quantitative verifica tion with 
experimental results. A sensitive method for measuring 
the residual stresses needs to be developed so that the 
effect of various grinding parameters could be evaluated. 
Further the theoretical model has been considered for 
single pass plunge-cut grinding. It needs to be extended 
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further to include the effect of successive grinning 
passes. The work could also he extended to include the 
effect of grinding wheel characteristics such as wheel 
hardness, grain size etc. 
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APPENDIX 


The elements considered in the solution of equation(3. 15 ) 
are sub parametric triangular elements as shown he low. 

j j Ov ^2.* 

s e /• \; _ ^ 

/ r" V' 6lC x ^) 

| ; «n- — \ 

6 : ^ .. 

where 1 , 2,3 are the primary node numbers and 4 , 5,6 are the 
secondary node numbers. 1 ^ , f 2 , ^3 are aatural coordinates 
of a point, Q, in the element. The variation of thermo elastic 
displacement potential, 4, within the element is taken as (38) 

i ''h N 

i < 

! * 2 ! ' 

1 i 

, 1 

, r 

; ^ ! 

. ^ e) = IN-, N 2 N 3 N 4 N 5 N 6 ] \ 5 { f (A.1) 

’ i 4 4 | 

| ^5 j 

i , j 

V_ 46 j> 

Here, the shape functions, til 1 , are given by 

N 1 = L 1 (2L 1 -1) ; N 2 = L 2 (2L 2 -1); ^3 = L 3 ( 21^-1) 

N 4 = 41^2 ; N 5 = 4 L 2 L 3 5 N 6 = ( ' A * 2 ' ) 

The Natural coordinates are related to the Cartesian coordinates 
(x,y) by the relationship (37) 
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■where are ' fc1ie Cartesian coordinates of node i and 
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the clement area, A v is given hy 
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The relations that are useful in simplifying the derivation 
of element stiffness matrix etc. are 
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The elements of the element stiffness matrix, [l4 e) ],and 
element load vector, l p (e) > . are given hy (38) 

T.T hi T>T 
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and p^ e) = -lid?. dA (3.19) 

Using equations (A. 2) through U-7) and equations (3.13) 
and (3.19), the elements of element stiffness matrix and 
element load vector can he derived as 


k (e) 
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Partial derivatives of the shape functions to he used in 
equations (3.23) are 
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The matrix of the derivatives of shapr functions to te used 
in equation (3.24) are 
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Table 1(a): Input data for calculating temperature distribution. 




•after Sauer (9). 

(Workpiece 

material 

nisi 1018 

Or steel 

Set 

No. 

v s 

ips 

ips 

d 

in. 

3? 

P 

lb£ 

h 

Btu 

m sec f 

a 1 

a 2 

a 

1 

770 

6.0 

0.00105 

13.4 

0.0000 

0.6912 

0.0000 

0.6912 

* 

770 

7.5 

O.OOI'I 

16.4 

0.0000 

0.4393 

0.1136 

0.5579 

3 

770 

7.95 

0.00056 

10.2 

0.0010 

0.7298 

0.0000 

0.7298 

4 

75 i 

9.5 

0.001 1 

29.9 

0.0000 

0.5143 

0.0086 

0.5229 

5 

770 

10.9 

0.0010 

30.5 

0.0000 

0.2172 

0.0974 

0.3146 

6 

750 

6.6 

0.0011 

13.8 

0.0000 

0.1420 

0.1366 

0.2786 

7 

750 

9.5 

0.00102 

21 .3 

0.0057 

a . 4 270 

0.0499 

0.4769 

8 

750 

12.1 

0.00099 

26.5 

0.0017 

0.3389 

0.0374 

0.3763 

9 

750 

9.2 

0.00098 

23.8 

0.0005 

0.3387 

0.0140 

0.3527 

10 

750 

9.2 

0.00102 

22.7 

0.0007 

0.4283 

0.0000 

0.4283 

11 

750 

9.6 

0.00095 

12.5 

0.0001 

0.6539 

0.0000 

0.6539 

12 

1360 

9.8 

0.00096 

7.0 

0.0001 

0.6200 

0.0000 

0.6200 

13 

750 

0.0270 

0.0305 

5.9 

0.0053 

0.1350 

0.0395 

0.1745 

14 

750 

0.0275 

0.0155 

2.59 

0.0028 

0.0073 

0.0873 

0.0951 

15 

750 

0.0570 

0.0310 

6.25 

0.0025 

0.0642 

0 . 0/- 27 

0.1069 

16 

750 

9.5 

0.00102 

21 .3 

0.0057 

0.4270 

0.0499 

0.4769 


Table 1(b): Results of the present work. 
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Table 2: The constants in equation (2.5). 
ter Sauer (9). 


Wheel 

type 

h 

rn 

c 

24H 

0.000007 

2.0 

60H 

0.0000147 

3.0 

BOH 

0.000029 

6.0 


Table 3: Ratio, r, of undeformed chip wi^th 
to chip thickness, -if ter Sauer (9). 


Grain size 

Ratio r 

36 

18 

46 

15 

60 

12.5 

SO 

10 

120 

7.6 
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FIG. 3 A TYPICAL TEMPERATURE DISTRI- 
BUTION IN SURFACE GRINDING. 
AFTER SAUER (9). 
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FIG. 15 A TYPICAL RESIDUAL STRESS DIS- 
TRIBUTION FOR GRINDING CONDITIONS 
SET No.1, TABLE (1) 
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